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Abstract 

The levels in 26 Na with single particle character have been observed for 
the first time using the d( 25 Na,p7) reaction at 5 MeV/nncleon. The measured 
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excitation energies and the deduced spectroscopic factors are in good overall 
agreement with ( 0 +l)hu; shell model calculations performed in a complete 
spsdfp basis and incorporating a reduction in the N — 20 gap. Notably, 
the 1 ^ 3/2 neutron configuration was found to play an enhanced role in the 
structure of the low-lying negative parity states in 26 Na, compared to the 
isotone 28 Ah Thus, the lowering of the lp 3 / 2 orbital relative to the O/ 7/2 
occuring in the neighbouring Z =10 and 12 nuclei - 25,27 Ne and 27,29 Mg - is 
seen also to occur at Z =11 and further strengthens the constraints on the 
modelling of the transition into the island of inversion. 

Keywords: 81V35, nuclear structure, 26 Na, 7 -Ray transitions, transfer 
reactions, shell migration 


The breakdown of the shell model magic number N — 20 for neutrons, 
in the case of neutron-rich isotopes near 31 Na and 32 Mg DEI. has led to 
the concept of an “island of inversion”, where neutrons preferentially occupy 
orbitals above the normal IV = 20 gap, leaving vacancies below it (3j 4], The 
mechanism is now understood in terms of the shell model to involve specific 
valence nucleon interactions HUE] as well as a monopole shift in the effective 
single particle energies [ 6 ]. In a wider picture, the monopole migration of 
levels is understood as arising from the tensor and three-body components of 
the nucleon-nucleon force between valence nucleons in partially filled orbitals 

BE]. 

The level schemes of nuclei with odd N, such as N = 15 and N — 17 
(Figure [I]), provide important insight into the evolution of the single parti¬ 
cle structure approaching the island of inversion. Between Z = 8 and 14, 
the proton d 5 / 2 orbital is filled, and this affects the energies of the neutron 
orbitals. The reference energy is the 7/2“ level, as it is expected that the 
energy of the / 7 / 2 neutron orbital (a j > orbital, i.e. j = £+ 1 / 2 ) will depend 
relatively weakly on the filling of the proton d 5 / 2 orbital (also j > [7]). In 
the case of 21,23 0 the J 77 = 7/2“ assignments are tentative [13] IS] but the 
small width observed in 23 0 [15] supports the 7/2“ assignment. The system- 
atics and also the shell model [15j suggest a near degeneracy between the 
3/2 + and 3/2“ states in oxygen. As may be seen, 26 Na falls in a key region 
where the shell closures at N = 20 and 16 are weakening and strengthening 
respectively. As such, compared to the isotone 28 Al, the levels in 26 Na corre¬ 
sponding to a neutron in the 0 ^ 3/2 orbital would be expected to lie higher in 
energy, whereas states with a neutron occupying the lp 3 / 2 orbital may be de- 


2 



generate with the O/ 7/2 levels or even lie below therrj^J I 11 order to explore this 
behaviour the present study has been undertaken to locate the correspond¬ 
ing levels using neutron transfer onto 25 Na, produced as a secondary beam. 
As will be shown, the detection of coincident gamma-rays has dramatically 
enhanced the excitation energy resolution and provided vital complementary 
information on the spins and parities of the populated levels. 

Previous studies have identified the ground state of 26 Na as having J n = 
3 + [T61HT] and there are three 1 + states fed by the beta-decay of 26 Ne [T8UT9] . 
Apart from the low-lying quartet of states (3 + , 1 + , 2 + , 2 + ) below 407 keV [ 20 ] . 
which have been very recently probed via Coulomb excitation [5T], the many 
observed states ; 2 Q. 221 23, US] have a structure that is almost completely 
unknown. The (d,p) reaction employed here will selectively populate the 
levels with a predominantly single-particle structure, which are those most 
useful for testing shell model predictions. 

The experiment was undertaken at the ISAC-II facility at TRIUMF. A 
pure beam of 5.0 MeV/nucleon 25 Na ions of intensity 3 x 10' pps was employed 
to bombard a self-supporting (CD 2 ) n deuterated polythene foil of thickness 
0.5 mg/cm 2 . The target was mounted at the centre of the SHARC silicon 
strip detector array m which was surrounded by 8 segmented germanium 
clover gamma-ray detectors of the TIGRESS array [25] . Four clovers were 
centred at 90° in the laboratory and four at 135°, giving a total absolute 
photopeak efficiency of some 3% at 1.33 MeV. In order to identify events 
in SHARC arising from fusion-evaporation reactions initiated on the carbon 
in the target (a source of significant background), an A1 foil followed by a 
thin plastic scintillator detector (the TRIFOIL [25]) were installed 400 mm 
downstream of the target. The A1 thickness of 30 //rri was chosen to stop 
fusion-evaporation residues whilst transmitting the beam and the products 
of direct reactions Further, the scintillator foil was thin enough (10 /jrn) 
to transmit the radioactive beam. 

Events were recorded whenever a particle was registered in SHARC, in¬ 
cluding a logic signal to indicate whether the TRIFOIL had fired during the 
same beam bunch. The TRIFOIL allowed the events corresponding to the 
production of 26 Na via the (d,p) reaction to be highlighted in the analysis 
in order to optimise the software gating. In the final determination of the 


3 A multiplet of levels for each neutron configuration is expected, according to the 
coupling of the neutron with the unpaired proton in the 0 ( 4/2 orbital. 
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absolute differential cross sections of the (d,p) reaction, the TRIFOIL re¬ 
quirement was not imposed, however, owing to its efficiency being dependent 
on the position on the foil, and hence on the angle of the recoil proton [28]. 
The coincident gamma-rays were recorded in TIGRESS and their energies 
were corrected for Doppler shift (v/c ~ 0.1) [281 . 

The excitation energy spectrum for states populated in 26 Na, recon¬ 
structed from the observed energy and angle of the protons, is shown in 
Figure [2] where the resolution is 350 keV (FWHM). This figure illustrates how 
closely spaced states could be distinguished using the gamma-ray data. The 
FWHM resolution in gamma-ray energy at 1.8 MeV, after Doppler correc¬ 
tion, was 18 keV at 135° and 23 keV at 90°. This resolution, rather than that 
of the proton-derived excitation energy spectrum, defined the precision with 
which individual states could be selected. On the other hand, the proton- 
derived excitation energy was critical in determining the energy at which the 
26 Na was populated in the (d,p) reaction (i.e. prior to gamma-decay of the 
states). I 11 addition, gating on this excitation energy was employed in the 
determination of the gamma-ray branching ratios [28]. 

Elastically scattered deuterons from the target were prominent in the 
spectrum of particle energy versus angle and the corresponding differential 
cross section was extracted [28] . Optical model parameters for d+ 26 Mg at 
an almost identical centre of mass energy [29] were employed to construct a 
theoretical angular distribution and this was found to reproduced well the 
form of the measured distribution (including the observed minimum). The 
best-fit normalisation thus gave the product of the deuteron target thickness 
and the total integrated beam flux. The uncertainty associated with this 
fitting was estimated to be 3%. An estimate of systematic uncertainties in 
this analysis was obtained by repeating the procedure with parameters taken 
from d+ 28 Si scattering at the same energy ffl. The variation was less than 
2 %. 

Figure [3] shows the Doppler corrected gamma-ray energy (A 7 ) plotted 
against the excitation energy in 26 Na (E x ) as derived from the proton en¬ 
ergy and angle. The superior resolution of the gamma-ray energy is clearly 
apparent, and the states that overlap in E x can be distinguished using 
the gamma-ray energy. A number of levels with ground-state gamma-ray 
branches (I? 7 = E x ) can be seen, including many of the more strongly pop¬ 
ulated states. The neutron separation energy is S n = 5.57 MeV [31] . The 
ground state branch generally has a lower and better defined underlying 
background than any other peak in Figure [3] and hence for the most part 
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only these peaks were chosen for initial analysi^ (Table 1). With suitable 
background subtraction [28], the yield of protons could be deduced for indi¬ 
vidual states in 26 Na as a function of laboratory angle. These distributions 
were converted to absolute differential cross sections by taking into account 
the geometry of SHARC, the elastic scattering normalisation, the measured 
gamma-ray branching ratios and the gamma-ray detection efficiency (cor¬ 
rected for Doppler and relativistic angular abberation effects). The overall 
systematic error associated with these effects was typically 5-6%. 

Differential cross sections for states populated in 26 Na are shown in Figure 
[4] plotted in terms of the laboratory scattering angle. These are compared 
with reaction calculations performed using the code TWOFNR [32] employ¬ 
ing the Adiabatic Distorted Wave Approximation (ADWA) of Johnson and 
Soper [ST] with standard input parameters [35] including the Chapel-Hill 
(CH89) nucleon-nucleus optical potential [36 ]. This formalism benefits in¬ 
ter alia from having no requirement for a deuteron-nucleus optical potential. 
The magnetic substate populations from TWOFNR were used to check that 
the gamma-ray coincidence requirement did not alter the shape of the dif¬ 
ferential cross sections by more than a few percent across the full range of 
proton angles [28] . Spectroscopic factors, S, were extracted by normalising 
the calculated differential cross sections to the data using the full range of 
angles shown in Figure |4} The transferred angular momentum, L, that best 
describes the shape of the measured angular distribution for each state is 
listed in Table 1. 

States for which more than one value is possible for the transferred an¬ 
gular momentum have been analysed by fitting a linear combination of the 
calculated cross sections using the two possible values Lq and L 2 (where 
L!=L). It is important to note that the contribution from the larger of hq 
and L 2 is suppressed (for a similar spectroscopic factor) due to poorer kine¬ 
matic matching. As a result, the spectroscopic factor deduced for the higher 
L is seen generally to exhibit a larger statistical uncertainty. As may be 
seen in Table 1, the addition of L 2 to the fit does not significantly change 
the spectroscopic factor deduced using just L 1; i.e. Si ~ S. The overall 
uncertainties assigned to the spectroscopic factors are some 20% (dominated 
overwhelmingly by the reaction calculations) [35]. The inferred spin and 


4 A full analysis including cascade decays and the more weakly populated states will 
follow. 


5 



parity assignments are discussed below. 

The 25 Na projectile has a structure, in the simplest shell model picture, of 
three 0d 5 /2 protons coupled to 5/2 + and a closed 0 d 5 / 2 neutron orbital. The 
positive parity orbitals lsi/ 2 and 0d 3 / 2 , as well as the negative parity orbitals 
O/ 7 / 2 , 1 ^ 3 / 2 ) lpi /2 • • • are thus available for neutron transfer. As shown in 
Figure [2j the known low-lying positive parity states below 450 keV [151 2D] 
are populated. As discussed below, the data also show clearly that there is 
population of negative parity states, as would be expected from the previous 
(d,p) studies of 25,27 Ne [12], 13], 27 Mg [37] and 28 A1 pH 07] {39]. In order to 
interpret the results, shell model calculations have been performed including 
all 0 hu sd-shcll configurations for positive parity states and all lhu excita¬ 
tions in an spsdpf basis for negative parity states. The program OXBASH 

[40] was used with the USD-A [32] and WBP-M [13] interactions respectively. 
The WBP-M interaction, which describes the neighbouring nuclei 25,2 'Ne and 
29 Mg in a consistent manner [13], is a modification of the WBP interaction 

[41] that shifts the energies of the /p-shell orbitals down by 0.7 MeV. It was 
previously noted [13] that for the neon and magnesium isotopes the strongly 
populated 3/2“ and 7/2“ states (S ~ 0.4-0.6) appear to track fairly closely 
the shift in the energy of the /p-shcll orbitals. In fact, earlier work by Bender 
and coworkers also applied a similar interaction called WBP-a that improved 
the agreement with the shell model for states in nearby (and only slightly 
higher Z ) isotopes of A1 [12] and P [43J ;H]. Excitation energies for the 
WBP-M and USD-A calculations were computed relative to the (positive 
parity) ground state calculated with the same interaction. The association 
of experimentally observed states with shell model states, indicated in Table 
1 and Figure [5] was based on the observed L-transfer, the excitation energy 
and the gamma-decay selectivity (Figure [6| as discussed below. 

Above ~2.5 MeV excitation, no states of positive parity are expected to 
be strongly populated, according to the shell model calculations. The states 
seen here at higher energies are therefore likely to have negative parity. Of 
these, the most easily identified is the strongly populated state at 3.511 MeV 
which is observed to have a yield dominated by L= 1 neutron transfer at 
large laboratory angles (Figure [4]). The only observed gamma-decay branch 
for this state is to the 3 + ground state, in line with the favoured decay pattern 
for the lowest 4“ state in the isotone 28 A1 [38]. This state is thus assigned to 
be the 4/ state in 26 Na. 

There are two strongly populated higher lying levels that gamma-decay 
via the 3.511 MeV level (Figure [6]). In terms of predicted levels in this 
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energy range, the obvious candidates are the lowest and states (Table 
1). These could reasonably be expected to gamma-decay via the 4“ state, 
according to the pattern observed [39 j in 28 Al. 

The decays of both the proposed and 6“ states proceed in part (20% 
and 31% respectively [25]) via the state at 2.116 MeV. Given that this level 
is populated via an L= 2 transfer (Figure [4]) and that it is fed by these higher 
lying negative parity states, it has a likely J n assignment of greater than 4, 
and a comparison of possible spins with shell model calculations indicates an 
assignment of 5 + (supported by both the energy and the weak spectroscopic 
factor). Indeed, there is an analogous state at 2.581 MeV in 28 Al, perhaps 
weakly populated in the decay of the 6^ state [32], that has a 5 + assignment 

USUIS]. 

The two reasonably strong states seen just below the 4~ are most natu¬ 
rally associated with the 2“ and 3“ states predicted in the shell model (see 
Table 1). In particular, the identification of the 3“ at 3.135 MeV is con¬ 
firmed by its gamma-decay branch to the lowest lying 1 + state at 1.507 MeV 
as occurs in 28 Al (Table 28.6 of ref. [46]). This is further supported by the 
mixed Li = 1 plus L 2 = 3 differential cross section (Figure [I]). 

The spectroscopic factors extracted from the proton differential cross sec¬ 
tions are compared to the shell model predictions in Table 1 and Figure [5] 
The lsi /2 strength that leads to 2 + states is concentrated in the two known 
states below 0.5 MeV and sums to approximately 0.45 in both theory and 
experiment. The tentative identifications of 3^ and 4 \ states are based on a 
comparison with the shell model calculations as follows. The observed 0c?3/2 
strength leading to J 71 " = 3 + is concentrated in one state near 1.8 MeV with 
a spectroscopic strength in very good agreement with the predictions. The 
0d 3 / 2 strength leading to the 4 + states has a comparable magnitude in the¬ 
ory and experiment and is located close to 2 MeV, so that the experiment 
is compatible with the strength being concentrated in the second 4 + state. 
We note that the USD-A interaction [32j predicts the excitation energies for 
positive parity states to within ~150 keV, in line with typical shell model 
accuracy. When we instead used the WBP-M interaction to calculate the 
energies of states with a strong 0d 3 / 2 neutron character (namely the 3 + and 
2 + states), then the predictions lay about 0.4 MeV lower than experiment. 
This is because the WBP-M incorporates the USD interaction m to com¬ 
pute the 0 hu positive parity levels and this is known to underestimate 0 d 3 / 2 
neutron energies [12]. 

The spectroscopic factors deduced for the 4~ state indicate a comparable 


Table 2: Mixing between lp 3 / 2 and 0 f 7 /2 configurations in 26 Na, compared with the isotone 
28 Al. Values of (2J+1)S from HQ are converted to S using the confirmed spins [46]. 




26 Na, present work 


28 Al, ref. [48] 


J w 

E 

Si(lp3/2) b) 

^2 ( 0 f 7 / 2 ) 

E 

Sl(lP3/2) S 2 (0f 7 / 2 ) 

3^ 

3.135 

0.06 0.01 

0.10 0.03 

3.591 

0.19 

0.60 

A- 

3.511 

0.25 0.01 

0.51 0.05 

3.465 

0.11 

0.62 

3^ 




4.691 

0.31 

0.08 

2 “ 




4.766 

0.41 

0.15 

2 ^ 




4.905 

0.24 

0.06 

3" 




5.134 

c) 


excitation energy in MeV 

b > statistical error 

in italics c -* beyond E x 

range of HE. 


strength to that predicted for the lp 3 / 2 transfer, albeit somewhat weaker. 
However, the shell model fails completely to reproduce the substantial O/ 7/2 
strength. The excitation energy is given accurately by the WBP-M cal¬ 
culation, as are those of the other negative parity states, within 100 keV 
(J 7r = 2 _ , 3 _ , 5 _ , 6 _ ). These latter states are all predicted to have a struc¬ 
ture that overlaps substantially with that of a neutron in the 0/ 7 / 2 orbital 
coupled to the 0d 5 / 2 proton of the 25 Na ground state. It is worthwhile noting 
that the 3“ state is predicted also to have a similar spectroscopic strength for 
the coupling with a lp 3 /2 neutron. Because of a better kinematic matching, 
the L= 1 transfer dominates the yield; a similar situation is observed for the 
lowest 3“ state in 28 Al as populated in (d,p) [SHI 09, 45]. In 28 A 1 , the 0/ 7 / 2 
measured spectroscopic factor actually exceeds that of the lp 3 / 2 orbital by 
a factor of three [48] (Table 2). Another 3~ state in 26 Na is predicted by 
the WBP-M calculations to have approximately equal mixing and to lie at 
4.462 MeV, whilst a predominantly lp 3 / 2 neutron state occurs at 4.774 MeV. 
These states may be populated but have not been identified in the present 
analysis. 

In the case of the 4~ state in 26 Na, the deduced lp 3 / 2 spectroscopic factor 
is twice as large as that reported for 28 Al whereas the 0/ 7 / 2 spectroscopic 
factor is effectively unchanged (Table 2). In contrast, for the 3~ state, the 
spectroscopic factors are three times and six times lower, respectively, than 
in 28 A1 [48, 46]. A clearer picture emerges from Table 2 if the relative magni¬ 
tudes of the IP 3/2 and 0 / 7 / 2 spectroscopic factors for each of the 26 Na states 
are compared with the behaviour in 28 Al. For both the 3“ and 4r states in 
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26 Na, the lp 3 / 2 spectroscopic factor is half the magnitude of that for 0/ 7 / 2 . In 
contrast, the analogous states in 28 A 1 exhibit lp 3 / 2 spectroscopic factors that 
are three to hve times smaller. This demonstrates an enhanced role emerging 
for the lpz /2 orbital in the structure of the low-lying negative parity states 
in 26 Na as compared to 28 Ah Indeed, in the case of 28 A1 the spectroscopic 
factors (Table 2) indicate that the states with predominantly lp 3 / 2 structure 
lie in the region of 4.8 MeV, significantly above the O/ 7 / 2 -dominated states 
which are closer to 4.0 MeV (Table 2 of ref. [48]). It would be very inter¬ 
esting if the spectroscopic factors for the higher lying 2 ” and states in 
26 Na could be measured and compared. As is evident in Figure [HJ a theory- 
based comparison between the lp 3 / 2 and Of 7/2 strengths in 26 Na indicates 
that the levels with a 1 p 3 / 2 structure are on average around 1 MeV below 
those with 0 / 7 / 2 structure (in fact, 0.83 MeV when weighted by the spectro¬ 
scopic factors), an inversion that is in accord with the systematics of Figure 

m 

In conclusion, as noted in the introduction and illustrated in Figure [lj 
the ordering of levels in nuclei with A ~ 25-30 evolves dramatically as they 
become more neutron-rich, driven largely by the interaction between protons 
in the Oc^s / 2 orbital and the valence neutrons. Importantly in this context, 
the results presented here confirm that the evolution is also manifest in 26 Na. 
In particular, the low-lying 3“ and 4“ states in 26 Na are found to exhibit an 
enhanced influence of the lp 3 / 2 neutron orbital, compared to the isotone 
28 Al. In addition the WBP-M shell model calculations, in which the fp- 
shcll orbitals are lowered to reduce the N = 20 gap by 0.7 MeV, succeed in 
reproducing the energies of the negative parity states in 26 Na as they do in 
the neighbouring neon and magnesium isotopes [HIS]. From a theoretical 
perspective, a less ad hoc description of the transition into the island of 
inversion represents an interesting and important challenge. 
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Relative level energy (MeV) 



Figure 1: (Colour online) Energies of lowest levels in isotones TV = 15 (dashed line) and 
TV = 17 (full line). Energies for TV = 13 are included for oxygen. The points for Z=10-16 
are derived from refs. HUM mi nn and compilations, whilst those for oxygen are 
from refs. [HI ITS] . The energy reference is chosen to be the 7/2 level (see text) and the 
Z position of 26 Na is indicated. 
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Figure 2: (Colour online) Excitation energy spectrum for states in 26 Na, reconstructed 
from the measured energy and angle of the protons from the d( 25 Na,p) reaction. The 
spectrum for all events is shown (red), plus the results for coincidences with the known 
gamma-rays of 233 keV (blue) and 407 keV (black histogram) in 26 Na. The data are for 
proton laboratory angles backward of 90°. A TRIFOIL requirement has been imposed in 
all cases (see text). 
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Figure 3: (Colour online) Plot of the excitation energy E x of 26 Na states populated via 
the d( 25 Na,p) reaction, as deduced from the proton energy and angle, versus the energy of 
any coincident gamma-ray (after Doppler correction), E 1 . Gamma-decays directly to the 
ground state lie along the diagonal E x = E 7 , whilst cascade decays result in events lying 
above the diagonal. Data are for laboratory proton angles backward of 90°. A TRIFOIL 
requirement has been imposed (see text). 
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Figure 4: Differential cross sections for the reaction d( 25 Na,p) at 5.0 MeV/nucleon. The 
results of ADWA reaction calculations are also shown, normalised to the data, for the 
angular momentum transfers indicated and listed in Table 1 (dashed line = fit with single 
angular momentum, labelled as L or as Li where another fit is also shown, full line = sum 
of two contributions Li and L 2 ). 
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Figure 5: (Colour online) Level scheme of 26 Na as deduced from the present work compared 
to the results of shell model calculations employing the USD-A and WBP-M interactions 
(see text). The lengths of the coloured lines correspond to the spectroscopic factors as 
follows: red s-wave, blue p-wave, green d-wave, orange /-wave. 
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Figure 6: (Colour online) Decay scheme of 26 Na as deduced from the present work, in¬ 
cluding the gamma-ray branching ratios. 
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